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COMPUTER SOLUTION OF UNSTEADY NAVIER-STOKES EQUATIONS 
FOR AN INFINITE HYDRODYNAMIC STEP BEARING 
by Henry A. Putre 
Lewis Research Center 

SUMMARY 

The two-dimensional, unsteady, incompressible, finite-difference Navier-Stokes 
equations, together with special inlet and outlet boundary conditions and corner bound- 
ary conditions, were programmed and solved for an infinite hydrodynamic step bear- 
ing. Where possible, the numerical bearing solution was compared to the Rayleigh an- 
alysis (one-dimensional) solution. The bearing length-to- clearance ratio for these cal- 
culations was fixed at 200. Calculations were carried out for Reynolds numbers (Re), 
based on wall velocity and upstream film thickness, ranging from 100 to 1000. 

The unsteady solution with impulse wall startup, which has not been previously pub- 
lished, was calculated for a Re = 100, optimum geometry (Upstream length/Downstream 
length = 2.45) bearing. This solution indicates that 99 percent of the steady-state hydro- 
dynamic lift is achieved in the time the wall moves one-half of the bearing length. In a 
practical bearing, this amounts to a lift startup time of about 0.0003 second, which is 
exceedingly small compared to practical bearing startup times. 

The steady-state solution was calculated for various Reynolds numbers and for an 
optimum and off-optimum (Upstream length = Downstream length) bearing geometry. In 
general, these steady-state calculations resolved the detailed velocity and pressure field 
structure in the step region. The important feature of the steady-state numerical solu- 
tion, which is not accounted for in the Rayleigh analysis, is the pressure drop due to 
fluid acceleration at the step (i. e. , inertia effect). This pressure drop was found to be 
as much as 47 percent of the peak pressure at the step at a Reynolds number of 1000. 

For the optimum geometry bearing, the lift from the numerical solution, including 
the effect of large step pressure drop, was virtually identical to the Rayleigh analysis 
prediction. For the off-optimum geometry, the step pressure drop in the numerical 
solution caused the lift to range from 2 to 14 percent below the Rayleigh analysis lift for 
Reynolds numbers from 100 to 1000. 

The numerical results are also presented as computer streamline plots and as a 
computer-generated fluid marker motion picture that is available on request. 



INTRODUCTION 


Recent advances in digital computers and in numerical solution techniques have in- 
spired the numerical solution of the full Navier-Stokes equations. The earliest numeri- 
cal solutions published were for steady incompressible flows and simple geometries such 
as the circular cylinder and the rectangular cavity (refs. 1 and 2). These were followed 
by the transient incompressible flow solutions, especially the vortex flow development by 
Fromm and Harlow (ref. 3). Vorticity and stream function were generally used as vari- 
ables, rather than velocity and pressure, for mathematical tractability. 

Subsequently, Harlow and Welch (ref. 4) published a transient incompressible solu- 
tion technique that solves directly for velocity and pressure. An important feature of 
this technique is that boundary conditions for complicated flows are treated in the prac- 
tical physical variables, and this technique could also be extended to three-dimensional 
problems. This transient technique was developed primarily for calculating free surface 
flows, but it can readily be applied to confined flows as well, as was done in a receipt 
paper by Donovan (ref. 5) for a square cavity. 

In this report the method, as used by Donovan, is further extended to a configuration 
with more complicated boundary conditions . The main objective of this work is to solve 
for the velocity and pressure distribution in an infinite hydrodynamic step bearing with 
no side leakage (see fig. 1) using the two-dimensional incompressible Navier-Stokes 
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Figure L - Model of infinite Rayleigh step bearing showing Rayleigh analysis pres- 
sure solution and actual pressure profile. 
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equations. This solution will require treatment of boundary conditions for long chan- 
nels with through-flow ends and a step. 

The infinite step bearing (also called the Rayleigh step) is of fundamental importance 
to the theory of hydrodynamic lubrication. In practical applications, such as journal and 
thrust bearings, the bearing is finite and some side leakage occurs. However, the infinite 
step is generally analyzed because it incorporates all other features of step bearings. 

The original analysis of this bearing was published by Lord Rayleigh (ref. 6) and was 
limited to a laminar one -dimensional treatment. This analysis neglected fluid acceler- 
ation (inertia effects) along the bearing and is strictly valid only for low Reynolds num- 
bers. Pressure profiles with and without inertia effects are sketched in figure 1. For 
the convenience of the reader, the Rayleigh analysis is given in appendix A. Constan- 
tinescu (ref. 7) and other investigators in a similar way have analyzed turbulent bearings. 
This report is concerned only with laminar flows to Reynolds numbers of 1000 or less. 
This is a practical operating range for gas bearings and in some cases also liquid bear- 
ings. 

The inertia effects in an inclined pad bearing were approximately calculated by 
Osterle and Saibel (ref. 8). They found the inertia effects to be a significant but small 
contribution to the bearing lift. For this same bearing, Lyman and Saibel (ref. 9) calcu- 
lated the lift during bearing startup. These analyses cannot be applied directly to the step 
bearing, and the laminar step bearing startup with inertia effects remains to be analyzed. 

In this report the two-dimensional results of a Navier-Stokes equation startup solu- 
tion for a Reynolds number of 100 and the steady-state solutions for Reynolds numbers 
ranging from 100 to 1000 are presented. These results give improved hydrodynamic lift 
estimates, including inertia effects, and a clearer understanding of the flow details in the 
infinite hydrodynamic step bearing. Important differences will be indicated between the 
pressure profiles of the one-dimensional Rayleigh solution and the present two- 
dimensional Navier-Stokes equation solution. 

A computer-generated fluid marker motion picture has been made for the Reynolds 
number of 100 bearing startup, and it is available on loan. A request card and a de- 
scription of the film are included at the back of this report. 


SYMBOLS 

h channel height (hg also called bearing clearance) 

L bearing length 

l length of computing mesh region 

p pressure above ambient in feed groove 
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dimensionless pressure, 




P h 2 

/iU(Lj +L 2 ) 


dimensionless pressure from Rayleigh analysis 
Reynolds number, Re = h^U /v 
dimensionless time, T = Ut/hj 
time step 
time 

bearing wall velocity in x-direction 
average x-direction fluid velocity in channel 
fluid velocity in x-direction 
fluid velocity in y-direction 


hydrodynamic lift per unit depth, W = 



p dx 


dimensionless hydrodynamic lift, 




W dimensionless hydrodynamic lift from Rayleigh analysis 

A#" increment in dimensionless lift 

x direction along bearing 

Ax, Ay space increments 
y direction across bearing 

p viscosity 

v kinematic viscosity, v = \ifp 

p density 

Subscripts: 

i mesh location corresponding to y 

j mesh location corresponding to x 

L left side value, found by extrapolating upstream values to the step 

R right side value, found by extrapolating downstream values to the step 

S at the step 
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1 upstream of step 

2 downstream of step 


ANALYSIS 

In this section the equations for incompressible flow in an infinite step bearing are 
given, and the bearing model shown in figure 1 is discussed. The third dimension (into 
the page in fig. 1), although important in many real applications, is assumed infinite and 
will be omitted in this analysis. In addition, important features of the computer codes 
are described. 


Navier-Stokes Equations 

For constant property incompressible flow of a Newtonian fluid, the two-dimensional 
equations of motion are 


5u | 3u J 

at 


auv 

ax ay 


1 dp + J 3 2 u 


dx 


\ 9 y 



(i) 


3v ( 3v + 3uv 

at ay ax 


lSp + Jd 2 v 


p 3y 


vax" 


A. 

ax ay j 


( 2 ) 


The continuity equation also applies and is 


iH = 0 

ax ay 


(3) 


Equations (1) and (2) are a rearranged form of the incompressible two-dimensional 
Navier-Stokes equations (p. 52 of ref. 10). 


Step Bearing Model 

Equations (1) to (3) can be solved in principle for the entire flow region, including 
the lubricant feed grooves as shown in figure 1. However, the typically long aspect ratios 
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Lj/hj and Lg/hg indicate that the flow over most of the bearing length outside the step 
vicinity will be one-dimensional. For this reason, and in order to conserve computer 
storage, the Navier-Stokes equation computing region extends over only a portion of the 
bearing length l ^ in the vicinity of the step (see fig. 1). The remainder of the flow 

outside the computing region is treated as one -dimensional. In agreement with standard 
practice, the bearing end pressures are assumed equal to the ambient feed groove pres- 
sure. (A more precise evaluation of the bearing end pressures, which is of secondary 
importance, would require a more extensive analysis including flow in the feed grooves.) 
The step bearing is therefore divided into three regions: 

One-dimensional unsteady flow with p(-L^) = 0 for 

-I^rsxrs-Zj 

Two-dimensional unsteady flow, Navier-Stokes equation solution for 

-Z 1 < x < l 2 

One- dimensional unsteady flow with p(+Lg) = 0 for 

1 2 ~ x ~ L 2 


The technique of choosing the computing lengths Z j and Z 2 is presented in a later 
section. 

The various step bearing configurations that were analyzed are now discussed. The 
configuration is specified completely by the four dimensionless quantities Re, 

(Lj + LgJ/hg, hj/hg, anc * L i/ L 2‘ The Re y n °l ds number has values ranging from 100 to 
1000. The aspect ratio, which in practice is fixed by fabrication tolerances for clear- 
ance and space availability for length, is fixed at (L^ + Lg)/!^ = 200 (typical values are 
hg = 0. 01 cm and + Lg = 2 cm). Two bearing geometries are analyzed. One, 
referred to as the optimum bearing, has proportions very close to the optimum Rayleigh 
bearing (see appendix A) with h^/hg determined by the numerical mesh size and 
Lj/Lg = (hj/hg)*^ according to equation (A6). The other is a typical off-optimum 
geometry. These bearing configurations are summarized in table I. The Rayleigh 
analysis step pressure and lift calculated from equations (A4) and (A5) are also shown. 
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TABLE I. - BEARING CONFIGURATIONS WITH RAYLEIGH 
ANALYSIS STEP PRESSURE AND LIFT 


Dimensionless 

quantity 

Optimum 

Off -optimum 

( L 1 + L 2^ h 2 

200 

200 

V h i 

11/20 

0.50 

V h i 

78. 1 

50 

V h i 

31.9 

50 

Re 

100, 300, 600, 1000 

100, 300, 600, 1000 

2? 

^s 

0.4120 

0.3333 

IT 

0.2060 

0. 1667 


Numerical Solution Technique 

Poisson equation for pressure . - The time -dependent Navier-Stokes equations and 
the continuity equation, as given in equations (1) to (3), are rewritten in a form amenable 
to numerical analysis. The basic technique used here is described in detail by Harlow 
and Welsh (ref. 4) and has been used by Donovan (ref. 5). Equations (2) and (3) are in 
turn differentiated with respect to x and y and then added to give 

, /~2 -.2 \ ,2 2 ,2 2 0 -,2 

l/3p i Sp\,3u , 3 v ,29 uv. 

For incompressible flow the bracketed term above vanishes and equation (4) becomes a 
time independent Poisson equation with spatial derivatives of velocities in the source 
term. The removal of time dependence with equation (4) is a key element of this nu- 
merical solution technique. 

The approach taken in the numerical solution is to consider the velocity values and 
derivatives in the source term of equation (4) known at the beginning of the n^ 1 time step 
and to solve the finite-difference Poisson equation for pressure at this time using suc- 
cessive over-relaxation. The pressure solution is thus calculated with a conventional 
numerical technique. With the pressure and velocity values known at the n** 1 time step 
after solving equation (4), the new velocities at the beginning of the (n + 1) time step are 
are calculated explicitly by using the finite difference forms of equations (1), (2), and (4), 
which are identical to those given by Donovan (ref. 5) and thus not repeated here. 


3 /9u + 3v^ 


= 0 


3t \3x 3y I 


(4) 


I 


Ax 



v l-(l/2),j 


Ay 

j u i, j-(l/2) 

' Pi,j ' 

u i,j+(l/2) 

T 

v iKl/2),j 



(a) Placement of field variables. 
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(b) Typical computing mesh layout 

Figure 2. - Typical field variable placement and computing mesh layout 


The placement of field variables about a mesh cell is shown in figure 2(a) . Pressure 
is only defined at the cell center and the u and v velocities are only defined on alternate 
cell sides. The finite difference mesh used for the Rayleigh step bearing calculation is 
shown in figure 2(b) which indicates that the top and bottom sides are solid walls and the 
left and right sides are through-flow ends. 

Boundary condition equations . - For a specific flow configuration, in addition to the 
general difference equations just described, realistic numerical boundary conditions must 
be carefully formulated compatible with the general equations. In this Rayleigh step 
bearing calculation the treatment at the open ends and at the step requires special atten- 
tion. Two types of boundary conditions need to be specified: 

(1) Initial values of u, v, and p at all points inside the mesh boundaries for the 

first time step 

(2) Boundary values of u, v, and p or certain spatial derivatives on the mesh 

boundary at every time step 

The initial values for the bearing startup are readily supplied since initially all u, 
v, and p are zero. For the steady-state calculations, a different Reynolds number 
steady-state solution could be reached while saving considerable computer time by start- 
ing with a previous solution and simply changing the viscosity. Thus the steady-state 
values from a previous calculation were used as initial values for a higher Reynolds num- 
ber calculation. 

The boundary values on the closed sides (B, D, E, and F of fig. 2(b)) are consid- 
ered next. No-slip velocity boundary values are used at these closed sides and the exter- 
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nal mesh values of u and v are chosen consistent with the terms in the momentum 
equation. Along the moving wall these u and v values are 


v ii,j = ' v 3§,j 

(5a) 

u 2,j+i =2U - u 3,j+! 

(5b) 

V 2i,j =0 

(5c) 


The pressure value outside the moving wall which is used in the Poisson equation 
solution is 


p 2,; 


= P 3 i (U ? ; .1 


Ax 






(5d) 


Note that equation (5d) follows directly from the momentum equation (eq. (2)) with the 
no -slip boundary conditions (eqs. (5a), (5b), and (5c)). If side B were stationary, these 
same equations with U = 0 in equation (5b) would apply. 

The boundary values on the open sides (A and C of fig. 2(b)) are complicated by the 
fact that the velocities in general cannot be independently specified as in equations (5a), 
(5b), and (5c), but are related to pressure through the time -dependent momentum equa- 
tions. Another important fact is that the open side boundary conditions be applied in such 
a way that the solution be least sensitive to the location of the boundary. Several formu- 
lations were used in trial computer rims. The boundary conditions on the in-flow side A, 
which resulted in least sensitivity of the pressure profile to the distance of side A from 
the step, are now given. 

The flow upstream and on the boundary is assumed one -dimensional so that the bound- 
ary conditions on side A are 


^P= — P (6a) 

9x Lj - l j 


= 0 (6b) 

dx 


dv 


0 


dx 


(6c) 
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Using equation (6b) in the continuity equation (eq. (3)) and integrating along side A 
give an additional restriction on side A: 


v = 0 (6d) 

When only equations (6a), (6b), and (6c) are used for numerical boundary conditions, 
the computer value of v on the boundary is generally nonzero, in disagreement with 
equation (6d). In order to make the momentum equation on side A exactly one- 
dimensional, equation (6d) was used as a numerical boundary condition rather than (6c). 
Equation (6c) will be used as an adjustment in the x-momentum equation. 

Equations (6b) and (6d) are written directly in terms of exterior mesh values 


u. ,1 = u. „l 

1 , J -2 1,02 


v i+i,2 " v i+|,3 


(7a) 

(7b) 


Equation (6c) cannot be applied in a similar manner since this would require a double 
valued v. a 0 . Instead this equation is applied directly in the numerical x-momentum 

1+2 j “ 

and Poisson equations. This is most easily done by adjusting equation (6a) with a term 
so that the derivative d v/dx dy vanishes in the x-momentum and Poisson equations. 
Thus the pressure boundary condition becomes 

dp , M a 2 v _ p 

3x dx dy Lj - 1 1 


The numerical boundary condition that combines equations (6a) and (6c) is 


P 


i,2 





7 i+i , 2 - v i- f , 3 + v i-i , 2^ 


(7c) 


The boundary values on the inflow side A are then calculated from equations 7(a), 
(7b), and (7c). A similar set of equations applies on the outflow side B. These open-end 
boundary conditions must be applied at sufficient distances from the step where the flow 
is very nearly developed. 


10 



I 


Corner Ambiguities 

At the 90° corners (as, e.g. , between sides A and B or E and F in fig. 2(b)), the 
velocity and pressure boundary values can be applied in a straightforward manner. The 
boundary conditions at the 270° corner (between sides D and E), however, are more dif- 
ficult to apply because some points on and outside the boundary are required to have more 
than one value. This situation was remedied by rewriting the finite -difference x and y 
momentum equations and the Poisson equation at the three cells that meet in this corner 
so that the boundary values are used directly in the equations. 


Computer Codes 

Numerical solution code. - The finite -difference equations and the boundary value 
equations were programmed in FORTRAN IV for calculations on an IBM 7094 computer. 

In the computer code the values of U, hj, and p are arbitrarily set at one. This code 
includes two accuracy checks which should be satisfied at each time step. The first, a 
convergence test, requires that the Poisson equation iterations give local pressure val- 

_ O 

ues that vary less than 10” /Re within successive pressure iterations. The second, a 

continuity and mathematical consistency check, requires that the quantity (9u/3x) + (3v/3y) 

_ o 

for each cell be small, less than about 2x10” /Re, if the incompressible difference equa- 
tions have been correctly applied. Time step values ranging from AT = 0.02 to 0.05 
were used in calculations. A typical mesh consisted of 100 mesh intervals of Ax = 0. 10 
by 20 mesh intervals of Ay = 0. 05. With a 100 by 20 mesh size and a time step of 
AT = 0. 05, about 90 minutes are required on the IBM 7094 to carry the calculations out 
to near steady state, or T = 46 at Re = 100. Listed output at specified time intervals 
from the code are u, v, p, continuity and pressure convergence checks at each cell, and 
integrated pressure along the moving wall. 

Visual display computer codes. - A second code was written for converting the large 
amount of information from the numerical solution code into a series of animated pictures 
which show the location of fluid markers as they are moved in time and space by the ve- 
locity field. The technique which results in a computed motion picture of the flow has 
been used by Harlow and Welsh (ref. 4) and also by Donovan (ref. 5). A motion picture 
was made for the Re = 100 startup case. Several frames from this movie are shown in 
figure 10 and will be discussed later (p. 19). 

A third code, which is basically similar to the movie code, was written to provide a 
streamline picture of the flow field at any time interval. This code dramatically shows 
in a single frame the very fine details of the flow field. Streamline pictures at very near 
steady state are shown for several cases in figure 9 (to be discussed later, p. 19). Of 
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particular note here is the fact that these pictures illustrate the flow pattern even in a 
region where velocities are very small, such as ahead of the step. 


RESULTS AND DISCUSSION 
Numerical Results 

The results of the numerical calculations are presented in several parts. First, a 
general accuracy check of this computer code is described. Next, the flow startup solution 
is discussed for the Re = 100 optimum bearing. This is followed by a discussion of the 
steady-state solution for Reynolds numbers ranging from 100 to 1000 for the optimum and 
off -optimum step bearings. Important differences between the numerical solution and the 
Rayleigh step solution, especially concerning pressure and hydrodynamic lift, are em- 
phasized. 

Comparison of code with previous work . - The computer code, described in the an- 
alysis, was temporarily modified with parabolic velocity profiles specified as the inlet 
and outlet boundary conditions and was used to calculate steady-state channel flow about 
an upstream facing step. The results were compared with an existing steady-state Navier- 
Stokes equations numerical solution at = 32 by Kawaguti (ref. 11). The velocity 

profiles calculated here and by Kawaguti are shown in figure 3. Good agreement was 
shown with reference 11 in the velocity profiles and in other important quantities such as 
pressure drop at the step. 

Computing region length for step bearing. - Preliminary computer runs were made 
with the boundary conditions given in the analysis section and with the open end boundar- 
ies at three distances from the step. A sufficient length was determined for the one- 
dimensional end boundary conditions by comparing solutions for short, medium, and long 
lengths ((z j/hj) + (Zg/hj) = 1+1, 2+2, and 3+7) for the Re = 100 steady-state solution. 
Compared to the long length solution, the medium length solution gave pressures that 
agreed within 0.6 percent, and the short length solution gave pressures within 1. 1 percent. 
The corresponding differences in lift were 0. 5 and 1. 1 percent. The end velocity pro- 
files were very nearly parabolic in these cases and were in closer agreement than the 
pressures. These results indicate that at Re = 100 the one -dimensional flow boundary 
conditions, applied at the distances z j = Z2 = 2hj from the step, give better than one per- 
cent pressure and lift agreement. The calculations, to be described, used the largest 
distances possible with a reasonable mesh size and the available computer storage, that 
is, z j = 3hj and 1 2 - 7h^. For the Re = 1000 calculation on a larger computing ma- 
chine, it would be desirable to lengthen the downstream computing region. 
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Figure 3. - Steady-state velocity profiles for channel flow with step calculation 
check for Re - (Uh/v) = 32. 


Unsteady solution with impulse wall startup, Re = 100. - The numerical solution for 
bearing flow startup was calculated for the optimum bearing at Re = 100. The velocities 
and pressures at every mesh point were calculated at time intervals of AT = 0.05 from 
rest to near steady state at T = 46 . The calculated inlet and outlet velocity profiles at 
various times are plotted in figure 4(a). Pressure distributions along the moving wall at 
various times are plotted in figure 4(b) . The time dependence of hydrodynamic lift during 
the flow startup is shown in figure 4(c) . 

The step bearing startup solution has not been previously published. Thus, these 
results, especially the lift-time dependence, contain new hydrodynamic lubrication infor- 
mation. For lack of an existing solution, the inlet and outlet velocity profiles are com- 
pared in figure 4(a) to the analytical Couette flow startup solution profiles (given in 
ref. 10, p. 74). Agreement at short times (T^4), when the step does not interfere 
with the moving fluid at the top wall, is excellent and constitutes a second accuracy 
check on the flow calculations. 

As shown in figures 4(b) and (c), the pressures on the moving bearing wall increase 
with time from startup. The pressure drop at the step in figure 4(b) is due to the fluid 
acceleration over the step. 
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Dimensionless pressure, 


Normalized 

time, 

T 



0 . 2 .4 .6 .8 , 1.0 0 .2 .4 .55 


Normalized distance from moving bearing wall, y/hj 

(a-1) Upstream of step, x/ hi = -3.0. (a-2) Downstream of step, 

x/hj = +7.0. 


(a) Velocity profiles at various times from bearing startup. 



Distance from step, x/hj 


(b) Pressure variation along moving wall during bearing startup. 

Figure 4. - Velocities, pressures, and lift calculated from impulse wall startup to steady state. 
Reynolds number, 100; L 2 = 2.45 Lj; optimum bearing configuration. 
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(c) Variation of hydrodynamic lift with time from bearing 
startup. 

Figure 4. - Concluded. 

The consequences of this step pressure drop are discussed in the next section. The 
lift variation with time shows that 99 percent of the steady-state value is reached in the di- 
mensionless time T = 46. In practical terms, this means that, with a moving wall speed 
of say 30 meters per second and a channel height hj of 0.02 centimeter, the bearing has 
developed 99 percent of its lift in a time of 0.0003 second. This flow relaxation time is 
very small when compared to practical wall startup times which are of the order of sec- 
onds. This indicates that the common practice, for unsteady wall speeds, of superim- 
posing steady-state lift solutions during any bearing transients is valid. 

Another interesting aspect of the lift-time curve is the sharp lift increase in fig- 
ure 4(c) at about T = 2. 5. This time delay to lift increase can be accounted for by using 
figure 4(a) as the time required for the upstream velocity to become nonzero at y = hg. 

A similar time delay was noted in reference 9 for the inclined pad bearing startup. 

Steady-state sol ution for various Reynolds numbers . - The same computer code that 
was used for the previous startup calculation was also used for steady- state calculations 
at various Reynolds numbers. The calculations for the optimum geometry are presented 
first, followed by the off-optimum geometry calculations. 

Steady-state velocity and pressure profiles were calculated for the = 2.45 Lg 
geometry at the Reynolds numbers of 100, 300, 600, and 1000. The calculations were 
carried out to A#/ AT < 10”^. At this cutoff, the velocity profiles on sides A and C are 
very nearly parabolic. Typical steady-state velocity profiles and the pressure along the 
moving wall are plotted in figures 5(a) and (b) for a Reynolds number of 300. The velocity 
profile plot in figure 5(a) indicates that beyond about one channel height upstream of the 
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Dimensionless velocity, u/U 


0 .2 .4 .6 .8 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 



(a) Typical velocity profiles at various distances from step. Reynolds number, Uh^/v = 300. 



Distance from step, x/hj Distance from step, x/h^ 


(b) Typical pressure variation along moving wall (Re = 300) over en- (c) Pressure variation along moving wall, in computing region, for 

tire bearing length. various Reynolds numbers. 

Figure 5. - Steady-state pressure and velxity profiles for optimum bearing for = 2.45 L 2 . 

step and three channel heights downstream of the step the velocity profiles undergo little 
or no change. This upstream distance was found to be nearly independent of Reynolds 
number, while the downstream distance was found to increase with Reynolds number. A 
narrow backflow region occurs near the stationary wall, upstream of the step, and is 
better illustrated in streamline plots that follow. 

The typical pressure variation is plotted in figure 5(b) for the entire bearing length. 
Both the numerical solution results and the Rayleigh analysis results are shown. The 
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numerical solution pressure drop in the vicinity of the step should be noted. This pres- 
sure drop is mainly due to fluid acceleration near the step. The net result of this step 
pressure drop is that in the upstream region the pressures and hydrodynamic lift are 
greater than predicted by the Rayleigh analysis. The opposite is true in the downstream 
region . 

A magnified plot of pressure at the moving wall, showing only the numerical calcu- 
lation region, is given in figure 5(c). Here the numerical results at various Reynolds 
numbers are compared to the Rayleigh analysis solution. From figure 5(c) it is seen that 
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(a) Typical velocity profiles at various distances from step. Reynolds number, Uh^/v 300. 




(b) Typical pressure variation along moving wall (Re = 300) over en- 
tire bearing length. 


(c) Pressure variation along moving wall, in computing region, for 
various Reynolds numbers. 


Figure 6. - Steady-state pressure and velocity profiles for off-optimum bearing for Lj = L 2 . 
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the step pressure drop effect increases with Reynolds number. 

The off-optimum geometry (Lj = Lg) calculations were carried out for Reynolds num- 
bers of 100, 300, 600, and 1000. Typical velocity profiles and pressure profiles are 
plotted in figure 6. A comparison of the velocity profiles in figures 5(a) and 6(a) shows 
that the Lj = L 2 geometry has a larger backflow region upstream of the step that oc- 
cupies about one-fourth of the channel height. This larger backflow region is also pre- 
dicted by the Rayleigh analysis. The moving wall pressures in the step region, shown in 
figure 6(b), have approximately the same variation as those for the optimum geometry in 
figure 5(b). 

An important difference between the numerical solution and the Rayleigh analysis has 
been pointed out: the pressure drop in the step region is not included in the Rayleigh an- 
alysis. This pressure drop will be evaluated by extrapolating the upstream and down- 
stream pressure values to the step to get the values of p^ and p^. This step pressure 
drop, normalized with p L , is plotted for both bearing geometries in figure 7 where it is 
seen to be nearly proportional to the Reynolds number. This dimensionless pressure 
drop varies from 7. 5 to 47. 3 percent for the optimum bearing and 5. 0 to 36.7 percent for 
the off-optimum bearing with Reynolds numbers from 100 to 1000. 

The steady-state hydrodynamic lift, which was calculated by numerical integration of 
the pressure at the moving wall over the entire bearing length, is plotted in figure 8 for 
the cases in figures 5(c) and 6(c). The lift is normalized so that the Rayleigh analysis 
values from table I appear in figure 8 as constants, independent of Reynolds number. The 



Figure 7. - Step pressure drop for various Reynolds 
numbers. (Bearing length)/Clearance = 200. 



Figured - Variation of hydrodynamic lift with Reynolds 
number. (Bearing lengthl/Clearance = 200. 
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Rayleigh analysis result that the optimum bearing has 24 percent more lift than the 
Lj = Lg bearing is shown. Comparing the Rayleigh analysis and the numerical solution 
results of figure 8, shows two important features. With the optimum bearing, in spite of 
sizeable pressure drops at the step which are not accounted for in the Rayleigh analysis, 
both solutions give identical lift values (within the numerical solution accuracy of about 
one percent). This surprising result is in agreement with the fact that, for the optimum- 
bearing only (where Ljhg = Lgh^ are given by eq. (A6)), the lift with an arbitrary abrupt 
pressure discontinuity at the step (given in eq. (A13)) is exactly the same as the Rayleigh 
analysis lift. The lift for the = L 2 geometry is seen in figure 7 to be significantly 
less than the Rayleigh analysis prediction. This reduction in lift is due to the pressure 
drop at the step and is approximately proportional to Reynolds number, with a lift reduc- 
tion of 1.9 to 14. 3 percent at Reynolds numbers of 100 to 1000. Therefore, depending on 
the bearing geometry, the step pressure drop may result in a significant correction, 
which increases with Reynolds number, that must be applied to the Rayleigh analysis lift 
calculation. 


Visual Flow Presentation 

As discussed in the analysis section, two types of computer flow visulization tech- 
niques were used with the numerical solution results to give a detailed picture of the flow. 
Both techniques utilized the numerical solution velocity values at one time interval to cal- 
culate the differential displacement of fluid marker particles during the time interval. 

The position of these markers after each displacement was plotted on a microfilm record- 
er display unit and photographed. 

Moving the markers with the steady-state velocity field and connecting the marker 
locations enables a single microfilm frame to present the complete streamline picture. 
These streamline pictures, with equal streamline spacing upstream, are shown in fig- 
ures 9(a) and (b) for the optimum bearing at Re = 100 and 300 and in figure 9(c) for the 
Lj = L 2 bearing geometry at Re = 100. Figures 9(a) and (b) show a vortex region for the 
optimum geometry that enlarges with Reynolds number. Figure 9(c) shows the larger 
size of the backflow region with the L 2 = geometry. Another interesting result in 
these streamline plots and in others for mesh spacings as small as 0. 025 by 0. 025 is that 
no vena- contra eta occurs along the top of the step. 

The other marker plotting technique which gives a new photograph after each set of 
displacements results in an animated motion picture that can be viewed in detail to show 
the time dependent nature of the flow. Figures 10(a) to (f) are a sequence of such motion 
picture frames during startup for the optimum geometry, Re = 100, case. 
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(a) Dimensionless time, T =0. 



(b) Dimensionless time, T = 5.0. 



(c) Dimensionless time, T = 10.0. 


Figure 10. - Computer motion picture frames for startup of optimum bearing. Reynolds number. 
Re = 100. (Upper wall moves from left to right) 




SUMMARY OF RESULTS 


The two-dimensional, unsteady, incompressible, finite -difference Navier-Stokes 
equations, together with special inlet and outlet boundary conditions and corner boundary 
conditions were programmed and solved for an infinite hydrodynamic step bearing, in the 
vicinity of the step. Where possible, the numerical bearing solution was compared to the 
Rayleigh analysis (one-dimensional) solution. In these calculations, the bearing length- 
to-clearance ratio (L^ + Lg)/!^ was fixed at 200. 

The important results of this study are as follows: 

1. A bearing startup solution for Re = 100 indicates that the step bearing with im- 
pulse wall startup achieves 99 percent of the steady-state lift in the dimensionless time 
Ut/hj = 46. This represents a practical startup time of the order of 0. 0003 second. 

2. The important feature of the steady-state numerical solution that does not appear 
in the usual Rayleigh analysis is the acceleration pressure drop at the step. This pres- 
sure drop for Reynolds numbers of 100 to 1000, expressed as a fraction of the peak pres- 
sure at the step, varies from 5. 0 to 36.7 percent for the off -optimum bearing and from 
7. 5 to 47.3 percent for the optimum bearing. 

3. The ultimate effect of the step pressure drop (inertia effect) on the hydrodynamic 
lift is less than the numerical accuracy of about 1 percent in the case of the optimum 
bearing. Therefore, the optimum bearing lift, but not pressure, can be accurately pre- 
dicted with the Rayleigh analysis. For the off-optimum bearing, the numerical solution 
resulted in values of lift that ranged from 1. 9 to 14. 3 percent below the Rayleigh analysis 
lift for Reynolds numbers from 100 to 1000. 

4. This report also presents the results of two computer flow visualization techniques. 
The first gives a series of streamline pictures which clearly show flow reversal ahead of 
the step. The second gives a series of animated motion picture frames that can be viewed 
in detail to show the time dependent nature of the flow . 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 26, 1969, 

122-29. 
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APPENDIX A 


RAYLEIGH ANALYSIS OF INFINITE STEP BEARING 

The equations generally used in the hydrodynamic lubrication theory are derived 
from these Navier-Stokes equations with certain simplifications. The simplifying as- 
sumptions in equations (1) to (3) that lead to the well known one-dimensional Reynolds 
equation of lubrication are listed in reference 12 (p. 6). The important ones for this 
analysis are the following: 

(1) There is no pressure variation across the film (Sp/3y = 0). 

(2) The flow is steady and fluid inertia is small compared to viscous shear (left 

sides of eqs. (1) and (2) vanish). 

2 2 

(3) The velocity term (3 u/3y ) dominates to the extent that the other viscous terms 
are negligible. 

Thus, the momentum equation becomes 


3X 3y 2 


(Al) 


This equation, together with the continuity equation, is integrated with no-slip velocity 
conditions at y = 0 and y = h and results in the one -dimensional Reynolds equation of 
lubrication for an infinite bearing with variable film thickness shape: 


_d 

dx 



= 6juU — 
dx 


(A2) 


This equation was solved by Lord Rayleigh (ref. 6) in a classic paper for several film 
shapes, and he showed that the stepped film gave the maximum lift of any film shape for 
fixed inlet and outlet film thicknesses hj and hg. 

The Rayleigh solution for the step bearing is straightforward and can be found in 
several texts, such as Pinkus and Sternlicht (ref. 12) and Cameron (ref. 13). It involves 
setting the right side of equation (A2) equal to zero for both halves of the step bearing and 
using the boundary conditions 


p(x = -Lj) = 0 

(A3a) 

p(x = 0") = p(x = O'*) 

(A3b) 

p(x = +Lg) = 0 

(A3c) 
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A typical Rayleigh solution pressure profile is shown in figure 1. The important 
Rayleigh analysis results of this solution are the peak pressure at the step 


p s = 


P B h 2 


M U(Lj +L 2 ) 


= 6 


L l L 2 h 2^ h l " h 2^ 


(Lj + L 2 )(L 2 h3 + Ljhg) 


(A4) 


and the hydrodynamic lift, where W = / 6 pdx 

• L 1 


/ +L 2 


ir = 


W / h 2 


= 3 


pUVLj +L 2y 


LiLgh^hi - h 2 ) 


(Lj + L 2 )(L 2 hj + Ljhg) 


(A 5) 


The bracketed quantity in equation (A5) is dimensionless and is seen to depend only 
on the geometric bearing proportions. A valuable use of this one -dimensional closed- 
form solution is that bearing geometry can easily be optimized to maximize the lift. For 
this optimization, the bearing length + L 2 and film clearance h 2 are considered 

fixed and the optimum ratios Lj/Lg and hj/h 2 can be found exactly by successive dif- 
ferentiation of equation (9) . 

The mathematical details are given in reference 12 (p.62), and the resulting equa- 
tions for the optimum ratios are 



V 


The optimum proportions that satisfy equations (A6) and (A7) are 

Vl + (AS) 

h 2 2 
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r=( 1+ T) V2 <A9) 

For this optimum bearing the dimensionless step pressure and the hydrodynamic lift of 
equations (A4) and (A 5) become 

p g = 0.4125 (A10) 

IT = 0.2063 (All) 

(Note that the difference between these values and those given in table I is due to the 
fact that the optimum bearing proportions in table I are adjusted to the nearest mesh 
spacing. ) 

An interesting result can be derived for the hydrodynamic lift in the case when a 
known pressure drop at the step is included in the analysis . The pressure continuity con- 
dition (eq. (A3b)) is replaced with the discontinuity 

p(x = 0“) = p L (A 12) 

p(x = 0+) = p R (A 13) 

When the Rayleigh analysis that led to equations (8) and (9) is extended, the results for 
an infinite step bearing with pressure discontinuity at the step are, for peak pressure at 
the step, 



and for hydrodynamic lift, 



(A 14) 


(A 15) 
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Note that, for the proportions given by equation (A6), the lift is the same as without the 
discontinuity at the step. 
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